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We briefly review the nonlinear sigma model approach for the subject of increasing interest: 
"two-step" phase transitions in the Gross-Neveu and the modified Nambu-Jona-Lasinio models 
at low N and condensation from pseudogap phase in strong-coupling superconductors. Recent 
success in describing "Bose-type" superconductors that possess two characterstic temperatures and 
a pseudogap above T c is the development approximately comparable with the BCS theory. One can 
expect that it should have influence on high-energy physics, similar to impact of the BCS theory on 
this subject. Although first generalizations of this concept to particle physics were made recently, 
these results were not systematized. In this review we summarize this development and discuss 
similarities and differences of the appearence of the pseudogap phase in superconductors and the 
Gross-Neveu and Nambu-Jona-Lasinio - like models. We discuss its possible relevance for chiral 
phase transition in QCD and color superconductors. This paper is organized in three parts: in the 
. first section we briefly review the separation of temperatures of pair formation and pair condensation 

in strong - coupling and low carrier density superconductors (i.e. the formation of the pseudogap 
phase ). Second part is a review of nonlinear sigma model approach to an analogous phenomenon in 
(f) ■ the Chiral Gross-Neveu model at small N. In the third section we discuss the modified Nambu-Jona- 

Lasinio model where the chiral phase transition is accompanied by a formation of a phase analogous 
to the pseudogap phase. 
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0^ ' Many concepts of particle physics have a close relation to superconductivity, for example the Nambu-Jona-Lasinio 
model §- @ was proposed in analogy to the BCS theory of superconductivity and is considered as a low-energy 
effective theory of QCD. Recently a substantial progress has been made in the theory of superconductivity in systems 
with strong attraction and low carrier density. That is, it has been observed that away from the limits of infinitcsimally 
weak coupling strength or very high carrier density, BCS-like mean-field theories are qualitatively wrong and these 
systems possess along with superconductive phase an additional phase where there exist Cooper pairs but no symmetry 
is broken due to phase fluctuations (the pseudogap phase). What may be regarded as an indication of the importance 
of this concept to particle physics is that recently the formation of the pseudogap phase due to dynamic quantum 
fluctuations at low N was found in the chiral Gross-Neveu model in 2 + e dimensions [Q . 
. ! Separation of the temperatures of the pair formation and of the onset of phase coherence (pair condensation) in 
strong-coupling superconductors, in fact, has been known already for many years (Crossover from BCS supercon- 
ductivity to Bose-Einstein Condensation (BEC) of tightly bound fermion pairs) |||(|. Intensive theoretical study of 
these phenomena in the recent years (see for example 0- f2l|| ), was sparked by experimental results on underdoped 
(low carrier density) cuprates that display "gap-like" feature above critical temperature T c that disappears only at a 
substantially higher temperature T*. There is experimental evidence that this phenomenon in high-T c superconduc- 
tors may be connected with precritical pairing fluctuations above T c . At present, this crossover has been studied by 
variety of methods and in many different models. 

Because of intimate relationship of many problems in particle physics to superconductivity it seems natural to guess 
that the pseudogap may become a fruitful concept in high energy physics too. 

Below we review these phenomena in superconductors and discuss its possible implications for QCD. 

II. PSEUDOGAP PHASE IN STRONG-COUPLING AND LOW CARRIER DENSITY THEORIES OF 

SUPERCONDUCTIVITY 
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I. INTRODUCTION 
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A. Perturbative results 



The BCS theory describes metallic superconductors perfectly. However, it failed to describe even qualitatively 
superconductivity in underdoped High-T c compounds. One of the most exotic properties of the latter materials is the 
existence of a pseudogap in the spectrum of the normal state well above critical temperature that from an experimental 
point of view, manifests itself as a significant suppression of low frequency spectral weight, thus being in contrast to 
the exactly zero spectral weight in the case of the superconductive gap. Moreover, spectroscopy experiments show 
that a superconductive gap evolves smoothly in magnitude and wave vector dependence to a pseudogap in normal 
state. Besides that, NMR and tunneling experiments indicate the existence of incoherent Cooper pairs well above 
T c . In principle it is easy to guess what is hidden behind these circumstances, and why BCS theory is incapable of 
describing it. Let us imagine for a moment that we are able to bind electrons in Cooper pairs infinitely tightly - 
obviously this implies that the characteristic temperature of thermal pair decomposition will also be infinitely high, 
but this does not imply that the long-range order will survive at infinitely high temperatures. As first observed in 
P|, long-range order will be destroyed in a similar way, as say, in superfluid 4 He, i.e., tightly bound Cooper pairs, at 
a certain temperature will acquire a nonzero momentum and thus we will have gas of tightly bound Cooper pairs but 
no macroscopic occupation of the zero momentum level q = and with it no long-range order. Thus phase diagram 
of a strong-coupling superconductor has three regions: 

• The superconductive phase where there are condensed fermion pairs. 

• The pseudogap phase where there exist fermion pairs but there is no condensate and thus there is no symmetry 
breakdown and no superconductivity. 

• The normal phase with thermally decomposed Cooper pairs. 

Of course, the existence of bound pairs above the critical temperature will result in deviations from Fermi-liquid 
behavior that make the pseudogap phase a very interesting object of study. In order to describe superconductivity in 
such a system the theory should incorporate pairs with nonzero momentum. Thus, the BCS scenario is invalid for 
description of spontaneous symmetry breakdown in a system with strong attractive interaction or low carrier density 
(see [^|, and references therein). So in principle in a strong-coupling superconductor onset of long range order has 
nothing to do with pair formation transition. The existence of the paired fermions is only necessary but not sufficient 
condition for symmetry breakdown. The BCS limit is a rather exotic case of infinitesimally weak coupling strength 
and high carrier density when the disappearance of superconductivity can approximately be described a as pair- 
breaking transition. The strong-coupling limit is another exotic case where the temperatures of pair decomposition 
and symmetry breakdown can be arbitrarily separated. There is nothing surprising in it: formally, in the case of Bosc 
condensation of 4 He we can also introduce a characteristic temperature of thermal decomposition of the 4 He atom; 
however this does not mean that this temperature is somehow related to the temperature of the Bose condensation 
of the gas of He atoms. A schematic phase diagram of a superconductor is in shown in Fig 1. 
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FIG. 1. Schematic phase diagram of a superconductor with arbitrary coupling strength. In the strong-coupling limit, 
the temperature of the superconductive phase transition tends to a plateau value corresponding to the temperature of Bose 
condensation of a gas of tightly bound fermion pairs, whereas the characteristic temperature of thermal pair decomposition 
grows monotonically as a function of the coupling strength. 

Let us show how one can obtain a pseudogap phase starting from the BCS Hamiltonian. This was first done in 
the pioneering work by Nozieres and Schmitt-Rink || and in the functional integral formalism for a system with 
5-function attraction by Sa de Melo, Randeria and Engelbrecht 0. In this subsection we briefly reproduce a part of 
the transparent article J^] and in the following section we will show how qualitatively the same result can be obtained 
within nonlinear sigma model (3D XY-model) approach proposed by the author Q. In the following sections we 
discuss analogous nonlinear-sigma model approach to the similar phenomena in the chiral GN and NJL models. 

The Hamiltonian of the BCS model is: 

F = E /d D ^(x)(-|^Wx)+ 5 |d D x^(x)^(x)^(x)^ T (x), (1) 

(7 

where Vv(x) is the Fermi field operator, a =t, I denotes the spin components, m is the fermionic mass, and g < the 
strength of an attractive potential gS(x — x). 

The mean-field equations for the gap parameter A and the chemical potential fj, can be obtained with a standard 
variation procedure: 

-- = -T J-tanh^, (2) 



n , . ....... 
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where the sum runs over all wave vectors k, N is the total number of fermions, V the volume of the system, and 

k 2 

2m 



E k = j£ + A 2 with & = (4) 



are the energies of single-particle excitations. 

The 5-function potential produces divergence and requires regularization. A BCS superconductor possesses a natural 
cutoff supplied by the Debye frequency u>d- For the crossover problem to be treated here this is no longer a useful 
quantity, since in the strong-coupling limit all fermions participate in the interaction, not only those in a thin shell 
of width ojd around the Fermi surface. To be applicable in this regime, we renormalize the gap equation in three 
dimensions with the help of the s-wave scattering length a s , for which the low-energy limit of the two-body scattering 
process gives an equally divergent expression S: 
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Eliminating g from (|^) and (^|) we obtain a renormalized gap equation 
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in which 1//cfo s plays the role of a dimensionless coupling constant which monotonically increases from — oo to oo as 
the bare coupling constant g runs from small (BCS limit) to large values (BEC limit). This equation is to be solved 
simultaneously with (^). These mean-field equations were analyzed e.g. in Ref. ||. 

In the BCS limit, the chemical potential n does not differ much from the Fermi energy 6f, whereas with increasing 
interaction strength, the distribution function broadens and a decreases, and in the BEC limit we have tightly 
bound pairs and nondegenerate fermions with a large negative chemical potential |/x| 3> T. Analyzing Eqns. (|^) 
and (H) we have from (||) the critical temperature in the BCS limit (/i ^> T c ) T C BCS — 8e~ 2 e~ > n~ 1 eF exp(— -k /2kF\a s \) 
where 7 = — r'(l)/r(l) = 0.577... , from (^) we have that chemical potential in this case is /1 = cf- In the 
strong coupling limit, from Eqs . ([§), (@) we obtain that in the BEC limit u = —Eb/2, where Ef, = 1/ma 2 is 
the binding energy of the bound pairs. In the BEC limit, the mean- field eq. (||) gives that the "gap" sets in at 
T* ~ Eb/2log(Eb/eF) 3 ^ 2 - A simple "chemical" equilibrium estimate (a b — 2/i/) yields for the temperature of pair 
dissociation: Tdi SSO c — Ebf log(£^/ei?) 3 / 2 which shows that at strong couplings T* is indeed related to pair formation 
H (which in the strong-coupling regime lies above the temperature of the onset of phase coherence ) . Obviously 
T* is a monotonous function of coupling strength. If we take into the account gaussian fluctuations we can see that 
in the strong coupling regime, the temperature T* , obtained with the above estimate is not related in any respect to 
the critical temperature of the onset of phase coherence. Expression for the thermodynamic potential with gaussian 
corrections reads |TJ||: 

n = fi -T^lnr(q,^). (?) 

q,i<3i 

where 

r-W) = T, l 1 -"^ + - (8) 



^ I iqi - £k - 6c+ q k 2 J Ana 



Where is Fermi occupation and iqi — iT2ttI. It is convenient following to ^ to rewrite fl in terms of a phase shift 
according to definition: r(q,w± iO) = |T(q, a;)| exp(±i#(q, w)). After inclusion of Gaussian correction the number 
equation N = —dfl/d[i reads: 



n = n {fi,T) + }^ — n B (u)) — (q,w) (9) 

q J — 00 "A* 



Where no is density of "free" fermions defined in (§) and ns(oj) = l/(exp(oj/T) — 1) is the Bose function. In order 
to study behavior of T c one should solve a set of the number and gap equations. In the BCS limit T c is not affected 
substantially by Gaussian corrections, thus the superconductive transition can be described by mean-field theory 
and correspondingly T c T*. Q In the opposite limit, numerical solution ||^| show that the temperature of the 
superconductive phase transition tends to a constant value that does not depend on the coupling strength and is equal 
to the temperature of condensation of the ideal Bose gas of particles of mass 2m and density n/2, where m and n are 
the mass and density of electrons correspondingly: 



2C(3/2) 



2/3 

— = 0.218e f (10) 
m 



Where 6f was used simply as a dimensional constant, namely the Fermi energy of the gas of free fermions with the 
density n and mass m (obviously at very strong coupling strength, when all fermions are paired there is no Fermi 
surface). 

The system of the gap and number equations can be solved analytically in the strong-coupling limit. First as it was 
pointed in one can make the following approximation: to retain Gaussian corrections only in the number equation 
and solve it together with mean-field "gap" equation. Near T c in the strong -coupling regime one finds that /x(T c ) = 
—Eb/2, where Eb is the energy required to break a pair. One can observe that in this limit T(q, z) has an isolated pole 



J As first discussed in 1960s, even in BCS superconductors there is a narrow region of precritical pairing fluctuations. This 
gives rise, e.g., to the so-called paraconductivity effect. In particle physics, this phenomenon was pointed out by Hatsuda and 
Kunihiro g|| 
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on the real axis for each q, representing a two-body bound state with momentum q. Since formally in this limit we can 
make energy required to break a pair arbitrarily large, this pole is widely separated from the branch cut representing 
the continuum of two-particle excitations. The low energy physics at temperatures much lower that temperature of 
thermal pair decomposition is governed by this pole and one can write T(q,iq m ) ~ R(c\)/[iq m — Wfc(q) + 2/i], where 
Wb(q) ~ — Eb + |q| 2 /4m. The partition function then may be written in the following form: 

z = z J d(j)d<j)exp <j ^2 <t>Sii - w &(q) + 2 m)0 9 > ■ (n) 

Correspondingly the strong-coupling number equation reads: 

n = n + n B [^h(q) - 2^]. (12) 
q 

From which neglecting n follows the result ( |l(i| ) |^.|| []. 



B. Nonperturbative nonlinear sigma-model (NLSM) approach to the BCS-BEC crossover in superconductors 

The above discussed crossover in the BCS model was recently studied in details perturbatively in a variety of 
approximations (see for example pO|JTl[ |). Qualitatively, essential features of this crossover can be reproduced in 
another simple model system - with the help of deriving of an effective non-linear sigma model (i.e. 3D XY-model) 
0. Moreover in the same framework of nonlinear sigma model one can study as well analogous crossover in 2D 
superconductor [p^| , pr|Jl9| fl that can not be addressed with discussed in the previous section perturbative method 
due to absence of the long-range order in 2D. In two dimensions T c is identified with a temperature of the Kosterltz- 
Thouless transition Tkt in 2D XY-model. In the strong-coupling (or low carrier density) regime Tkt lies significantly 
lower than the temperature of the pair formation [p^| , pT|Jl9| j7| . 

As we have shown in many essential features known from numerical study of strong-coupling and low carrier 
density superconductors are reproduced with very good accuracy within 3Z?AY-model approach. In particular, in 
the framework of NLSM approach there is no artificial maximum in T c in the regime of intermedium couplings (which 
appears in the approximation discussed in the previous section). 



1. XY-model approach to 3D superconductors 

Let us now reproduce results of previous subsection in the framework of NLSM approach proposed by the author 
§• 

It is very transparent to study the properties of the BCS-BEC crossover to employ "modulus-phase" variables. 
Following to Witten p4j we can write a partition function as a functional integral over modulus and phase of the 
Hubbard-Stratonovich field Aexp(iy>): 

Z(fi,T) = Jl>AV(pexp[-pSl(T,A(x),d<p(x))]. (13) 

Assuming that phase fluctuations do not affect local modulus of the complex Hubbard -Stratonovich field we can write 
thermodynamic potential as a sum of "potential" and "kinetic" (gradient) terms: 

n(A(x), d<p(x)) ~ n grad (T, a, d<p(x)) + n pot (T, a) = / d 3 .? J(T ' A) (v y ) 2 + n pot (r, a). (u) 



2 The above estimate by |]|)| when we retained corrections only in the number equation renders correct limiting result ( |ic| ) 
however as it was observed by the authors of in the intermediate coupling re gim e it gives an artificial maximum in T c as 
a function of coupling strength, this artifact is removed in higher approximations |lfj 
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Obviously in the above expression the effective potential n pot (T, A) coincides with the ordinary mean-field effective 
potential. The gradient term fl gra d(T, A, dip(x)) coincides with the Hamiltonian of 3D XY model with a stiffness 
J(A,T). 

Let us reproduce low-temperature expression for the phase stiffness in strong-coupling regime from : 



J = ^L_^y™ T 3/2 exp 

4m 167T 2 



T 



(15) 



Where n and m are density and mass of fermions. We can see that thermal corrections to the first term in this regime 
are exponentially suppressed and the r.h.s. tends in this limit quickly to 

n 

Jbe = -—- (16) 
4m 

The form of this expression is not surprising - we see that at sufficiently strong coupling strength all fermions are 
bound to pairs and stiffness becomes equal of the low-temperature phase stiffness of the Bose gas of density n/2 
and boson mass 2m. Obviously all information about internal structure of composite bosons is evaporated from this 
expression in this approximation since at low temperature in this regime there are no thermal pair decomposition 
effects []. Thus we see that low-temperature expression for the stiffness of the phase fluctuations reaches a plateau 
value with increasing coupling strength, whereas temperature of the thermal pair decomposition is a monotonously 
growing function of the coupling strength. 

In principle knowledge of lowest gradient term governing gaussian fluctuations is not sufficient for the study of the 
position of the phase decoherence transition in a system with preformed pairs. In continuum, 3D XY model is a free 
field theory and there is no phase transition. Phase transition of 3D XY model was studied in great details on a 
lattice, so we can consider a lattice model of BCS-BEC crossover and we can verify aposteriory to what extend lattice 
model reproduces features of this crossover in continual model. Many aspects of relation of 3D XY model to Bose 
condensation, in particular derivation of the Gross-Pitaevskii equation near T c can be found in [ p6[ . 

So, let us consider theory ([IJ) on a lattice with spacing a = l/n^ir where n pair is concentration of Cooper pairs. 
This model describes condensation of hard-core bosons on the lattice and is a special case of Bose-Hubbard model ^. 

Critical temperature of the phase transition of 3D XY model can be obtained with a simple mean-field estimate 

U& 

Tl DXY « 3Ja (17) 

— 1/3 

Where a = n pa / r is the lattice spacing. In contrast to the ordinary 3-PXy-model, in order to find temperature of the 



phase transition we should study a system of the equations (^), (||) , (15) and (p~7[) . System of these equations can be 
solved however analytically in the strong coupling limit, the result is 



\2J n V3 2V6^3/2 J c m exp 



2m 

With increasing coupling strength, this quickly tends from below to the value (compare with (|i"o|)): 



(18) 



3n 2 / 3 3 

^w^-wrr***"*' 2 **- (19) 

Where constant ep is Fermi energy of free fermi gas of density n and fermion mass m. We should observe that in 
nonperturbative NLSM approach T c approaches plateau value (J19|), that depends only on mass 2m and density n/2 



3 In Q we employed a finite-temperature generalization of the gradient expansion at T = discussed in B| . 

4 In the ordinary attractive Hubbard model considered in jf| the critical temperature is a nonmonotonous function of the 
coupling strength: T c decreases in the strong-coupling limit due to in that model composite bosons move via virtual ionization. 

We should stress that estimation of critical temperature of the phase transition of the effective 3D XY model with mean-field 
methods has nothing to do with BCS mean field approximation since derivation of the phase stiffness coefficient ( Jl5| ) required 
studying the gaussian fluctuations in the BCS model Q, and thus this may be regarded as the approximation of the same level 
as considered in the previous section. 
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of composite bosons, from below. This is in agreement with numerical study in higher approximations JU^, whereas 
in the approach presented in previous subsection T c has an artificial maximum at the intermediate coupling strength 
thus approaching the limiting value from above. Another circumstance that we discuss below is that NLSM approach 
gives also qualitatively correct results in the opposite limit of weak-coupling strength Q . In the weak-coupling limit 
near T c , the stiffness coefficient may be derived with the help of Gorkov's well-known method : 



This is precisely the coefficient of the gradient term in the Ginzburg-Landau expansion. In the weak-coupling limit 
the two temperatures of the onset of pairing correlations and the onset of phase coherence T* and T c merge according 
to the formula 0: 

T C = T* — i— J ^ -> T*, (21) 

With it one can see |?j that in the weak-coupling limit, the temperature of the phase transition of the effective 
XY-model tends from below to the characteristic temperature of the disappearance of the effective potential and 
merges with it for infinitcsimally weak coupling strength. Thus we arrive at some sort of aposteriori verification of 
BCS behavior in this limit in the model of hard-core composite bosons on the lattice (i.e. in this limit if nonzero 
modulus of the complex gap function Ae liptyX ^ appears at some temperature, at the same temperature phase coherence is 
established and continuous symmetry is broken) . In the weak and moderate coupling strength limits the disappearance 
of superconductivity is a competition of two processes - pairbreaking which is thermal excitation of individual particles 
and decoherence process which is thermal excitation of collective modes. 

Let us now summarize the results that follow from the NLSM consideration. In this model in strong-coupling or 
low carrier density regimes system possesses three phases: 

1. Superconductive phase (T < T 3DXY ). 

2. Pseudogap phase (T 3DXY < T < T*) - the phase where there exists a local gap modulus A that signalizes 
existence of the tightly bound (but noncondenced) fermion pairs but phase is random so average of the complex 
gap is zero (< |A| exp(i0) >= 0). So in this phase there is no superconductive gap and with it no symmetry 
breakdown. 

3. Normal phase (T > T*) the phase with thermally decomposed Cooper pairs. 



2. XY-model approach to 2D superconductors 

In two dimension there is no proper long-range order and superconductive phase transition is associated with a 
Kosterlitz-Thouless transition. In order to study this transition it is sufficient to extract lowest gradient term that 
determines temperature of the phase transition according to the formula |l3) , |l4| []l9| , fj 

Tkt - \ J(M, T kt , A( M , Tkt)). (22) 

This equation just like in the discussed above 3D case should be solved self-consistently with equations for the gap 
modulus and chemical potential (||), (^). The result in the strong-coupling limit is 

7r n f 1 

Thus for increasing coupling strength, the phase-decoherence temperature Tkt tends very quickly towards a constant 
value |]J|,[h|,[l9],[7j corresponding to KT transition in system of bosons with density n/2 and mass 2m (whereas 



(23) 



In principle there is no KT phase transition in a charged system due to Meissner effect, however coupling to electromagnetic 
field is always neglected in discussion of 2D superconductors, due to experimentally in-plane penetration length in high-T c 
materials is much larger than coherence length. 
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characteristic temperature of the thermal pair decomposition continues to grow monotonously with the growing 
coupling strength): 

Tkt = \~- (24) 
8 m 

So, this phenomenon in 2D is qualitatively similar to the above discussed 3D case. 



There is no superconductivity in the pseudogap phase however it exhibits rich exotic non-Fermi-liquid behavior 
due to local pairing correlations that makes it as interesting an object of theoretical and experimental study as the 
superconductive phase itself. In particular, along with specific heat, optical conductivity and tunneling experiments 
there are following circumstances observed in the pseudogap phase: In experiments on YBCO a significant suppression 
of in-plane conductivity o~ a b(u) was observed at frequencies below 500 cm -1 beginning at temperatures much above T c . 
Experiments on underdoped samples revealed deviations from the linear resistivity law. In particular cr a b(u> = 0;T) 
increases slightly with decreasing T below a certain temperature. NMR and neutrons observations show that below 
temperatures T* much higher than T c , spin susceptibility starts decreasing. 

In conclusion, let us once more emphasis essential features of this phenomenon in superconductivity: 

• Away from a very special limit of infinitesimally weak coupling strength and high carrier density, superconductors 
are characterized by two temperatures T c and T*(>> T c ). T* is the characteristic temperature below which 
pair correlations become important (or in the regime of strong interaction it is the characteristic temperature 
of the formation of real bound pairs). T c corresponds to the onset of phase coherence in a system of preformed 
fermion pairs. The region of non-Fermi liquid behavior between T c and T* calls the pseudogap phase, however 
the term "pseudogap" , originated in early experimental papers, may seem somewhat misleading since, even 
though a substantial depletion of low-frequency spectral weight is observed in this region experimentally - there 
is no superconductive gap in the spectrum. 

• One should note that there is no proper phase transition at T*, which is simply a characteristic temperature 
of thermal decomposition of certain fraction of noncondensed Cooper pairs. Even though the position of this 
temperature may be reasonable estimated with mean-field methods, second-order phase transition at T* is 
certainly an artifact of the discussed above approximation. Experiments on specific heat indicate however 
certain features at this characteristic temperature. 

In what follows we discuss possible implication of these results to QCD that may posses a phase analogous to the 
pseudogap phase in strong-coupling superconductors. The simplest model related to particle physics that displays 
pseudogap behavior of dynamical origin is Chiral Gross-Neveu model at low N that is discussed in the next section 
in 2 + e dimensions. 



III. PSEUDOGAP PHASE IN CHIRAL GROSS-NEVEU MODEL IN 2 + e DIMENSIONS AT LOW N 

Let us now discuss a phenomenon similar to pseudogap in a simple field-theoretic model - chiral version of the 
Gross-Neveu model |27j], whose Lagrange density is 

C = $ a i$lp a + ^ {^a^a) 2 + {^ailb^a) 2 • (25) 

Where index a runs from 1 to N. Appearance of the pseudogap phase in this model has quite similar roots with the 
above discussed phenomenon in strong-coupling superconductors. In superconductors we observed appearance of the 
pseudogap phase on the phase diagram in the region away from the limits of infinitesimally weak coupling strength 
or extremely high carrier density - i.e. in the regime when BCS mean-field treatment is no longer valid. Chiral 
Gross-Neveu model can be treated in the limit of infinite number of field components TV in a framework of mean-field 
approach quite similar to BCS theory. Transparently in the mean-field approximation one can find only one phase 
transition at certain value of the coupling strength, similar to BCS phase transition. However, at low N, system start 
to perform dynamic chiral fluctuations which as we have shown in give rise to a second, phase disorder, transition. 
So at low N the model possesses two transitions at two characteristic values of renormalized coupling constant. Let 
us reproduce this result. 
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One can write the collective field action for this model as 

f 



A co ii[cr} — N ^ — 7T— (c 2 + 7r 2 ) — iTrlog [i$ - a{x) - i l5 ir] 
z 9o 



(26) 



This expression is invariant under the continuous set of chiral 0(2) transformations which rotate a and ir fields into 
each other. This model is equivalent to another one: 



Here C is the matrix of charge conjugation which is defined by 



C'fC 



"7 



(27) 



(28) 



In two dimensions, we choose the 7-matrices as 7 = a , 7 1 = — ia 2 , and C — 7 . The second model goes over 
into the first by replacing ip — > — (1 — 75)^ + — (1 + r y5)Crp T , where superscript T denotes transposition. In the 



Lagrange density ( p7| ) we introduce a complex collective field by adding a term [N/2go) 
the partition function 



N 



, leading to 



Z[rj, fj]= J VtpV^VAVA^ exp < 1 



d u x 



Ipaifilpa + I (AtyjfCtya + C.C.) + IpT) + fj ij) - ^-\A\ 2 

2 250 



(29) 



The relation with the previous collective fields a and tt is A = a + in. 

Following to BCS procedure we can fix phase of the order parameter, then for a constant A, the effective action 
gives rise to an effective potential that in 2 + e dimensions reads: 



A 2 

9oV e 



2+e 



where fj, is an arbitrary mass scale, and the constant b e stands for 



b e = ^ 2 S D T(D/2)T(l - D/2) = l_L_.r(l - D/2), 



(30) 



(31) 



which has an e-expansion b e 
introduced by the equation 



so that 



— [l — (e/2) log (2ne 7 )] /ne + 0(e). A renormalized coupling constant g may be 
1 
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1 



1-1* 



AT v ' 2 

Extremizing this we obtain either Ao = or a nonzero Ao that solves the gap equation 

d D p 1 



l=. 9o tr(l) J 



(2n) 2 p 2 + A 2 ' 



in the following form: 



l _9^^£ fA 
2 



9 



(32) 



(33) 



(34) 



(35) 



where g* = —l/b e « 7re. 

In the limit N — > 00 this result is exact. On the other hand in the opposite limit of low N the system starts to 
perform fluctuations around saddle-point solution and in order to describe this system properly one should go beyond 
mean-field approximation and study propagator of the 6*-field - where 6 is the phase of the order parameter. 
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Let us consider first the case e = where the collective field theory consists of complex field A with 0(2)-symmetry 
A = |A|e iei . Such a system possesses macroscopic excitations of the form of vortices and antivortices that attract each 
other by a logarithmic Coulomb potential. It is known that in a such field theory involving a pure phase field 
9{x), with a Lagrange density 



(36) 



where j3 is the stiffness of the 9- fluctuations, there is a Kostrelitz-Thouless transition when the stiffness falls below 
/3kt - 2/tt |2|. 

Let us return to Gross-Neveu model. Performing an expansion around saddle point solution one can find propagator 
of the 0-field when e = MM: 



<-*ee ~ T7^7 + regular terms. 
TV q z 



Comparing this with the propagator for the model Lagrange density (|3 

r 1 1 



(37) 



(38) 



we identify the stiffness (3 = N/Air. 

The pair version of the chiral Gross-Neveu model has therefore a vortex-antivortex pair breaking transition if TV 
falls below the critical value TV C = 8 [f|. 

Consider now the model in 2 + e dimensions where pairs form when renormalized coupling constant becomes larger 
than the critical value g — g* « ire. In this case the expression for the stiffness of phase fluctuations reads 0] 



What implies a KT transition in the neighborhood of two dimensions Q at: 

9* 



N r . 



811-2- 
9 



(39) 



(40) 



The resulting phase diagram is shown in Fig. 
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FIG. 2. The two transition lines in the TV — g-plane of the chiral Gross-Neveu model in 2 + e dimensions. In order to stress 
difference between the local gap (i.e. "pseudogap") and the order parameter analogous to the "superconductive" gap we denote 
by M modulus of the order parameter (M = |Ao|). In this model M plays the role of the "quark" mass. For e = 0, the vertical 
transition line coincides with the TV-axis, and the solid hyperbola degenerates into a horizontal line at TV C = 8. In the limit 
TV — > oo generation of the quark mass happens simultaneously with "phase ordering" transition. 



7 There is a misleading statement about 3D case in 
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In the chiral formulation of the same model in 2 + e dimensions, the "pseudogap" phase has chiral symmetry in spite 
of a nonzero spontaneously generated "quark mass" M = |Ao| ^ 0. This phase is directly related to the pseudogap 
phase of a strong-coupling superconductor - where there are Cooper pairs but there is no symmetry breakdown due 
to violent phase fluctuations. The reason why this is possible is that the "quark mass" depends only on |Ao|, thus 
allowing for arbitrary phase fluctuations preserving chiral symmetry. It is very easy to sec that the solid hyperbola 
in Fig. ^ is not simply the proper (albeit approximate) continuation of the vertical line for smaller N. There are two 
simple arguments. One is formal: For infinitesimal e the first transition lies precisely at g = g* = ire for all N, so that 
the horizontal transition line is clearly distinguished from it (stiffness of the phase fluctuation in the regime <?*/<?—> 0, 
just like in the case of superconductor, reaches plateau value that does not depend on the coupling strength and e). 
The other argument is physical and also has a clear analogy in the corresponding phenomena in superconductivity. If 
N is lowered at some very large g, the binding energy of the pairs increases with 1/N {in two dimensions, the binding 
energy is 4Msin 2 [7r/2(A r — 1)]}. It is then impossible that the phase fluctuations on the horizontal branch of the 
transition line, which are low-energy excitations, unbind the strongly bound pairs. Accuracy of the "BCS" scenario 
in the limit N — > oo is clearly seen from the form of phase stiffness which has a factor N that "freezes" the phase 
fluctuations in this limit and thus all the physics is essentially governed by the size of the gap modulus. 

The 2 + 1-dimensional Chiral Gross-Neveu model [^9| also exhibits an analogous behavior at finite temperature 
where a similar effect is governed by thermal fluctuations. At finite N the temperature of KT transition deviates 
from mean-field temperature of the gap modulus formation, however in D=2+l the phase diagram is substantially 
different from the phase diagram of the same model inD = 2 + eatT = (see detailed discussion in ^] )[]. 



IV. CHIRAL FLUCTUATIONS IN THE THE NJL MODEL AT ZERO TEMPERATURE 



31 1 to generalize to the NJL model the nonlinear-sigma approach for description 
The authors 



Recently an attempt was made 
of chiral fluctuations proposed in 4jj7j] . The authors ||31[ claimed that at N c — 3 the NJL model does not display 
spontaneous symmetry breakdown due to chiral fluctuations. We show below that the NLSM approach does not allow 
to prove that chiral symmetry is always restored by fluctuations in the NJL model at N c — 3. Below we also discuss 
differences from the chiral GN model, where the NLSM approach allows one to reach a similar conclusion at low N. 

The Lagrangian of the NJL model reads M 



C = ipifiip - 



go 
2N n 



(lptp) 2 + (ipXai^Y 



(41) 



The three 2 x 2-dimensional matrices A a /2, generate the fundamental representation of flavor SU(2), and are normal- 
ized by tr(A Af,) = 25 a b- One can introduce Hubbard - Stratonovich fields a and 7r a : 



C = ijj (i$ - g - i7 5 A a 7r a ) ij) - — ^ (cr 2 + 7T 2 ) 



(42) 



After integrating out quark fields, following a standard mean- field variation procedure one can choose the pseudoscalar 
solution 7r a to vanish and the scalar solution a = M to be given by a gap equation: 



1 ./ ,w f d °P 1 

- = z(tr / l)(tr 7 l ) y {2iT)Dp2 _ M2 



(43) 



The momentum integral is regularized by means of a cutoff A. The constituent quark mass M in the limit N c — > oo 
is analogous to the superconductive gap in the BCS limit of the theory of superconductivity. 

At finite N c one can study fluctuations around the saddle point solution. The quadratic terms of expansion around 
the saddle point are: 



It should be emphasized that the existence of the pseudogap phase in 2+1 dimensions due to thermal fluctuations in the 
Chiral Gross-Neveu model can not be rigorously proven in contrast to 2 + e-dimensional case discussed in this section (when there 
are two small parameters e and 1/N). The 2+1 dimensional problem lacks a small parameter which would allow to estimate 
accurately the position of T* at very small N, however it can be argued that rough low-iV estimates show the appearance of 
thermal-fluctuations-induced pseudogap phase in this model, especially pronounced at N < 4. 
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a'{q) 



G- 




Go 



a'(-q) 



(44) 



where (a' , V Q ) = (cr — Af, 7r a ) and G CT ^ are the inverse bosonic propagators. Implementing a momentum cutoff A, we 
can write G~\ for small qs as: 



g: 



"(27T) 



In 1 



A 2 



A 2 



M 2 / A 2 + M 2 



-Z(M/K)q%; 



g: 



-Z(M/A)(q 2 E + 4M 2 ) 



(45) 



In analogy to 31? AY"-model approach to strong-coupling superconductivity the authors of |3lJ] introduced a unit 
vector field rii = (no,n a ) = (a,n a )/p and set up an effective nonlinear sigma-model 



a r l P 
A K = 7T 



d 4 x[an,(x)] 2 . 



(46) 



The prefactor /3 = M Z(M/K), that follows from Eqs. (|4J) and (45) plays the role of the stiffness of the unit field 
fluctuations. 

Now let us observe that from the arguments given in |n| it docs not follow that the NJL model necessarily remains 
in a chirally symmetric phase at N c — 3. At first, in contrast to the (2 + e)-dimensional case, discussed in [Q, one 
can not make unfortunately, any similar calculations in a closed form in 3 + 1-dimensions because this theory is not 
renormalizable. It was already observed in |3^]- (34| that the cutoff of meson loops cannot be set equal to the cutoff 
for quark loops and thus the 1/N C corrected theory Q possesses two independent parameters that may be adjusted 
at will. We present another argument of a different nature rooted in the nonuniversality of the critical stiffness of a 
NLSM in four dimensions, which does not allow one to reach the conclusion of |n| in the framework of the NLSM 
approach. Our observation also applies to the NLSM description of precritical fluctuations in general systems. It also 
allows us to show that the additional cutoff discussed below can not be related to the inverse coherence length of the 
radial fluctuations in the effective potential as suggested in (3^j3l|] . 

The authors of by deriving G Wi(T have essentially extracted two characteristics from the initial system: the 
stiffness of the phase fluctuations in the degenerate minimum of the effective potential and the mass of the radial fluc- 
tuation. However, knowledge of these characteristics does not allow one in principle to judge if directional fluctuations 
will destroy long range order or the system will possess a BCS-like phase transition. The reason is that the critical 
stiffness of the nonlinear sigma model is not an universal quantity in 3 + 1-dimensions. So in principle knowledge of 
the stiffness of NJL model is not sufficient for finding the position of the phase transition in the effective nonlinear 
sigma model. The situation is just like that in a Heisenberg magnet, where the critical temperature depends on the 
stiffness along with lattice spacing and lattice structure. Thus if one is given only a stiffness coefficient one can not 
determine the temperature of the phase transition. The situation is in contrast to the 2D case where the position of 
a KT-transition can be deduced from the stiffness coefficient |28| . 

Let us recall a procedure for expressing the critical stiffness of the 0(4)-nonlinear sigma model via an additional 
parameter: one can relax the constraint n 2 — 1 and introduce an extra integration over the lagrange multiplier A, 

rewriting Eq. ©as: (f3/2) J d 4 x {[drii(x)] 2 + A [n 2 (x) - l] }. Integrating out the n, (^-fields, yields: 



A [A] =-f3 J d A x^- + ^Trln [-d 2 + X(x) 



(47) 



where N n is the number of components of rii(x) and Tr denotes the functional trace. This yields a gap equation: 

, d 4 k 1 

= N n 



(2tt) 4 k 2 + \' 



(48) 



The model has a phase transition at a critical stiffness that depends on an unspecified additional cutoff parameter 
that should be applied to the gap equation: 



/3 cr - N„ 



d 4 k 1 
(2tt) 4 fc 2 ' 



(49) 



For example, in the case of magnets the additional cutoff needed in Eq. ( |49| ) is naturally related to the lattice spacing. 
In p5|] a criterion was proposed that states that one can relate the inverse coherence length extracted from radial 
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fluctuations in an effective potential of an initial theory to the cutoff in the integral (^) so that all the parameters 
in a theory would be expressed from quantities derived from an initial model, and thus this modified model possesses 
a universal critical stiffness. However, unfortunately there is no reason for relating the cutoff needed in Eq. ( |49| ) 
to the coherence length of the modulus fluctuations and moreover we show that this procedure leads in general to 
unphysical consequences. It should also be observed that it does not make the theory consistent anyway because 
of the following circumstance: The authors of |u| suggested that, "since pions in the symmetry broken phase are 
composite they are not "defined" over the length scales much shorter than the inverse binding energy of the pair wave 
function which is equal to 2M" . Thus the authors of [|TJ performed integral in ( [49| ) up to cutoff AM 2 . However, since, 
as suggested in |3l[], the pion fields are not "defined" over the length scales much shorter than the inverse binding 
energy this may only serve as an estimate for "upper boundary" of what would be the universal critical stiffness value. 
So, unfortunately one can not make the conclusion of absence of symmetry breakdown in such a modified theory by 
observing that stiffness derived from the initial model is smaller than the maximal possible value of would be universal 
critical stiffness. It was also supposed in that the relation of the coherence length to the cutoff in the equation 
( [49l ) yields a universal criterion for judging the nature of symmetry breakdown in general physical systems. There is 
a simple counterexample: in the case of a strong-coupling superconductor, the effective nonlinear sigma model that 
describes fluctuations in a degenerate valley of the effective potential is a 3D XY-model. In the continuous case it is 
a free field theory and has no phase transition at all. The phase transition appears only in the lattice theory and, of 
course its temperature depends on the lattice spacing. As we discussed above with increasing coupling strength the 
low-temperature phase stiffness of the effective 3D XY model tends to a plateau value J = n/Am, where n and m are 
the density and the mass of fermions Jjj . Thus the temperature of the phase transition of the effective 3D XY-model 
is 

T 3DXY K n 
m 

where a is the lattice spacing. To be careful one should remark that accurate analysis shows that a strong coupling 
superconductor possesses two characteristic length scales: size of the Cooper pairs which tends to zero with increasing 
coupling strength, and a coherence length that tends to infinity with increasing coupling strength as the system evolves 
towards a weakly nonideal gas of true composite bosons |^J2^] . First if one relates the constant a in ( p0[ ) to the size 
of the Cooper pairs following to the arguments of j3lj] one will come to an incorrect conclusion of the absence of 
superconductivity in strong-coupling superconductors, in a similar way as the authors of |MJ] came to the conclusion 
of the nonexistence of symmetry breakdown in the NJL model. This is in a direct contradiction with behavior of the 
strong coupling superconductors discussed above. Second, if one attempts to relate a in ( |50|) to the second length 
scale of the theory, namely, the true coherence length, which tends to infinity with increasing coupling strength, then 
one will also come to a qualitatively incorrect conclusion |pq| . 

Thus, in general, the nonlinear sigma model approach for precritical fluctuations possesses an additional fitting 
parameter, which is the cutoff in the gap equation (fl9|), which can not be related to the inverse coherence length 
extracted from radial fluctuations in an effective potential. Thus within the NLSM approach one can not prove if the 
NJL model displays necessarily the directional fluctuations-driven restoration of chiral symmetry at low N c . 



V. CHIRAL FLUCTUATIONS AT FINITE TEMPERATURE AND A MODIFIED NJL MODEL WITH A 

PSEUDOGAP 

This section is based on the paper (3^]. The authors of |3l| employed NLSM arguments in an attempt to show 
that the NJL model cannot serve the study of the chiral symmetry breakdown. We have shown in the above that this 
conclusion appears to be incorrect since the critical stiffness in 3+1-dimensions is not a universal quantity and one 
has an additional fitting parameter. This is an inherent feature of the discussed NLSM approach in 3+1 dimensions 
(compare with the cutoffs discussions in nonrenormalizable models in a different approach |Q, and also |3^)). The 
above circumstance allows one to fix the critical stiffness from phenomenological considerations. However, we argue 
below that, what is missed in is that, in principle, the low-iV c fluctuation instabilities, when properly treated, have 
a clear physical meaning. Moreover, we argue that one can employ a NLSM for describing the chiral fluctuations (e.g. 
at finite temperature), provided that special care is taken of the additional cutoff parameter. Indeed, it was already 
discussed in the literature that at finite temperatures the chiral phase transition should be accompanied by developed 
fluctuations ( H||| and references therein). We argue that this process at low N c should give rise to a phase analogous 
to the pseudogap phase that may be conveniently described within a nonlinear sigma model approach. There are 
indeed other ways to describe these phenomena, however the NLSM approach seems to be especially convenient. The 
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description of the two-step chiral phase transition and appearance of the intermediate phase requires one to study the 
system at the next-to-mean-field level. Unfortunately, the NJL model is not renormalizable and does not allow one 
to make any conclusions about the importance of fluctuations in a closed form Q] . On the other hand, a pseudogap 
phase is a general feature of Fermi systems with composite bosons. The NLSM construction discussed below, because 
of its nonperturbative nature, can not be regarded as a regular approximation but may be considered as a tractable 
modification of the NJL model that has a pseudogap. One can also find an additional motivation for employing these 
arguments in the fact that the NLSM approach allows one to prove the existence of the phase analogous to pseudogap 
phase in the chiral GN model ||,[30) which is the closest relative of the NJL model. Also the NLSM approach works 
well for the description of precritical fluctuations in superconductors [Q - where essentially the same results have been 
obtained with different methods and in different models. We stress that these phenomena are a general feature of any 
Fermi system with attraction. Also, to a certain extent similar crossovers are known in a large variety of condensed 
matter systems. In particular, besides superconductors we might mention the exitonic condensate in semiconductors, 
Josephson junction arrays, itinerant and local-momentum theories of magnetism and ferroelectrics. 

Let us now consider the chiral fluctuations in the NJL model at finite temperature. Then, following standard 
dimensional reduction arguments [see e.g. f39fl l, the chiral fluctuations should be described by a 3D 0(4)-sigma 
model. Thus one has the following gap equation for the effective NLSM (i.e. the finite temperature analog of (fig)): 

J T f d 3 k 1 

The temperature of the phase transition of the three dimensional classical 0(4) sigma model with stiffness Jt is 
expressed via the additional parameter At needed in (|5l| ) as : 

2 A T 

The stiffness of thermal fluctuations Jt can be readily extracted from the NJL model. At finite temperature the 
inverse bosonic propagator of the collective field it for small q can be written as: 



G- 1 - -2 D ' 2 N [ V 
* C J (2tt) 3 ^ 



T 



(p 2 + M 2 +oj 2 l ) 



q 2 



2 D / 2 N r 



d 3 p 



1 1 /#TF\ 1 1 , - a ( yV + M 2 

8 (p 2 + M 2 ) 3 / 2 \ 2T J 167V + M 2COS I 2T 



q 2 



K(T, A T , M, N c )q 2 , (53) 



where A^ is a momentum cutoff. The propagator (|53|) gives the gradient term that allows one to set up an effective 
classical 3D 0(4)-nonlinear sigma model: 

J T (T,A T ,M,N C ) f , 3 r ^ 



E = 1K 1 y ' L; J d'xldn^x)] 2 , (54) 
where 

J T (T, A T , M, N c ) = K(T, A T , M, N c ) M 2 (T, A T ) (55) 

is the stiffness of the thermal fluctuations in the degenerate valley of the effective potential. The temperature - 
dependent quark mass M that enters this expression is given by a standard mean-field gap equation which also should 
be regularized with the cutoff Ay: 

I = 2x2 c / 2 V / f fP 2 J 2 (56) 

.90 (27r) 3 p 2 + M 2 + to 2 

It can be easily seen that, when we approach the temperature T* where the mass M{T) becomes zero, the stiff- 
ness J(T, At, M, N c ) also tends to zero. Formula d54| ) defines a generalized Heisenberg model with a temperature- 
dependent stiffness coefficient. The position of the phase disorder transition in a such system should be determined 
self-consistently by solving the system of equations for T c and M(T C ). Apparently just as in a superconductor with a 
pseudogap, the phase transition in such a system is a competition between the thermal depletion of the gap modulus 
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(this roughly corresponds to thermal pair breaking in a superconductor) and the process of thermal excitations of the 
directional fluctuations in the degenerate minimum of the effective potential. The "BCS" limit corresponds to the 
situation where T* merges with T c and it is easily seen that this scenario always holds true at N c — > oo. That is, at 
infinite N the mean-field theory is always accurate just as BCS theory works well in weak-coupling superconductors. 
In the framework of this NLSM construction, at low N c the scenario of the phase transition depends on the choice of 
M(0), At, and At, which should be fixed from phenomenological considerations. 



VI. CONCLUSION 



The precursor pairing fluctuations is a general feature of any Fermi system with composite bosons and it is the 
dominant region of a phase diagram of strong-coupling and low-carrier-density superconductors. At the moment it is 
a subject of increasing interest in different branches of physics. In the first part of this paper we briefly outlined the 
nonlinear sigma model approach to this phenomenon in superconductors in two and three dimensions. In the second 
part we discussed similar phomena in the chiral Gross-Neveu and Nambu-Jona-Lasinio models. This discussion should 
have relevance for hot QCD and color superconductors. 

We also note that in some sense similar phenomena are known in a large variety of condensed matter systems, in 
particular, besides superconductors we can mention itinerant and local-momentum theories of magnetism, exitonic 
condensate in semiconductors, ferroelectrics and Josephson junction arrays [ ^0[ . 

We would like to stress that the main purpose of this paper is to summarize present discussions of precursor fluctu- 
ations in the Gross-Neveu and Nambu-Jona-Lasinio models. We illustrated the discussion with a few examples from 
superconductivity, outlining occurrence of similar phenomena in several arbitrarily chosen models of superconductors 
with pseudogaps. Thus this paper can not be regarded as a review of this phenomenon in superconductivity which 
evolved recently to a very large branch of condensed matter physics. Thus our references to the papers on supercon- 
ductivity are by definition incomplete, for more complete set of references reader may consult corresponding reviews 
on superconductivity [e.g. the review |P3| ]. 
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